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Abstract 

Possible phase transition of strongly interacting matter from hadron to a Quark-Gluon Plasma 
(QGP) state have in the past received considerable interest. It has been suggested that this 
problem might be treated by percolation theory. The Color String Percolation Model (CSPM) 
is used to determine the equation of state (EOS) of the QGP produced in central Au-Au collisions 
at RHIC energies. The bulk thermodynamic quantities- energy density, entropy density and the 
sound velocity- are obtained in the framework of CSPM. It is shown that the results are in 
excellent agreement with the recent lattice QCD calculations(LQCD). 
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1. Introduction 

One of the main goal of the study of relativistic heavy ion collisions is to study the 
deconfined matter, known as Quark-Gluon Plasma(QGP), which is expected to form at 
large densities. It has been suggested that the transition from hadronic to QGP state 
can be treated by percolation theory [1] . The formulation of percolation problem is con- 
cerned with elementary geometrical objects placed on a random d-dimensional lattice. 
The objects have a well defined connectivity radius A, and two objects can communicate 
if the distance between them is less than A. Several objects can form a cluster of com- 
munication. At certain density of the objects a infinite cluster appears which spans the 
entire system. This is defined by the dimensionless percolation parameter ^ [2]. Percola- 
tion theory has been applied to several areas ranging from clustering in spin system to 
the formation of galaxies. Figure 1 shows the transition from disconnected to connected 
system at high densities. 

In nuclear collisions there is indeed, as a function of parton density, a sudden onset 
of large scale color connection. There is a critical density at which the elemental objects 
form one large cluster, loosing their independent existence. Percolation would correspond 
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Fig. 1. Left panel: Disconnected discs. Right panel: Overlapping discs forming a cluster of communication. 

to the onset of color deconfinement and it may be a prerequisite for any subsequent QGP 
formation. 

The determination of the EOS of hot, strongly interacting matter is one of the main 
challenges of strong interaction physics (HotQCD Collaboration) [3]. Recent LQCD cal- 
culations for the bulk thermodynamic observables, e.g. pressure, energy density, entropy 
density and for the sound velocity have been reported [3]. A percolation model coupled 
with hydrodynamics has been utilized to calculate these quantities from the STAR data 
at RHIC energies in central Au+Au collisions. 

2. Color String Percolation Model 

Multiparticle production at high energies is currently described in terms of color strings 
stretched between the projectile and target. Hadronizing these strings produce the ob- 
served hadrons. The strings act as color sources of particles through the creation of qq 
pairs from the sea. The number of strings grows with the energy and with the number 
of nucleons of participating nuclei. Color strings may be viewed as small discs in the 
transverse space filled with the color field created by colliding partons. Particles are pro- 
duced by the Schwinger mechanisms [4]. With growing energy and size of the colliding 
nuclei the number of strings grow and start to overlap to form clusters [5,6]. At a critical 
density a macroscopic cluster appears that marks the percolation phase transition. 2D 
percolation is a non-thermal second order phase transition, but in CSPM the Schwinger 
barrier penetration mechanism for particle production and the fluctuations in the asso- 
ciated string tension due to the strong string interactions make it possible to define a 
temperature. Consequently the particle spectrum is "born" with a thermal distribution 
[7,8]. The percolation threshold at which the spanning cluster appear, a "connected" 
system of color sources, identifies the percolation phase transition. 

With an increasing number of strings there is a progression from isolated individual 
strings to clusters and then to a large cluster which suddenly spans the area. In two 
dimensional percolation theory the relevant quantity is the dimensionless percolation 
density parameter given by [5,6] 
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where N is the number of strings formed in the colhsions and is the transverse area 
of the single string and S'at is the transverse nuclear overlap area. The critical cluster 
which spans Sn, appears for > 1-2 [9]. As ^ increases the fraction of Sn covered by 
this spanning cluster increases. 

The percolation theory governs the geometrical pattern of string clustering. It requires 
some dynamics to describe the interaction of several overlapping strings. It is assumed 
that a cluster behaves as a single string with a higher color field corresponding to the 
vectorial sum of the color charge of each individual string. Knowing the color charge, 
one can calculate the multiplicity /i„ and the mean transverse momentum (pf of the 
particles produced by a cluster of strings. One finds [5,6] 
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(Pt)n = ^^{p',)l (3) 

where ni and {pf)i are the mean multiplicity and average transverse momentum squared 
of particles produced by a single string with a transverse area Si = Trrg. In the saturation 
limit, all the strings overlap into a single cluster that approximately occupies the whole 
interaction area, one gets the following universal scaling law 

{Pt)n = -^ —IJ-n (4) 

This scaling law shows a reasonable agreement with all the experimental data for all 
projectiles, targets and energies [10] . In the limit of high density one obtains 

= VF\0 (5) 

with 



no = y^^ (6) 

being the color suppression factor. It follows that 

M = ArF(0m,(p2) = _L_(p2)^ (7) 

A similar scaling is found in the Color Glass Condensate approach (CGC)[11]. The satu- 
ration scale Qs in CGC corresponds to {Pt)i/F{^) in CSPM. The net effect due to -F(^) is 
the reduction in hadron multiplicity and increase in the average transverse momentum of 
particles. The CSPM model calculation for hadron multiplicities and momentum spectra 
was found to be in excellent agreement with experiment [12,13]. Within the framework of 
clustering of color sources, the elliptic flow, V2, and the dependence of the nuclear mod- 
ification factor on the azimuthal angle show reasonable agreement with the RHIC data 
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Fig. 2. Percolation density parameter f vs Npart 

[14] . The critical density of percolation is related to the effective critical temperature and 
thus percolation may be the way to achieve deconfinement in the heavy ion collisions [7] . 
An additional important check of this interacting string approach was provided by the 
measurement of Long Range forward backward multiplicity Correlations (LRC) by the 
STAR group at RHIC [15]. 



3. Experimental Determination of the Percolation density Parameter ^ 

To obtain the value of ^ from data, a parameterization of p-p events at 200 GeV is 
used to compute the pt distribution [16] 

dN,/dp^t ^a/ipo+ptr (8) 

where a, po, and a are parameters used to fit the data. This parameterization also can be 
used for nucleus- nucleus collisions to take into account the interactions of the strings [16] 

f {nSi/Sn)Au-AuY^'^ ,„^ 

Po -> Po / Q /Q \ (9) 

V {nSi/bn)pp J 

where S'„ corresponds to the area occupied by the n overlapping strings. The thermody- 
namic limit, i.e. letting n and 5„ — >■ oo while keeping ^ fixed, is used to evaluate F(^) 
given by Eq.(5). For nucleus-nucleus coUisions Eq.(8) becomes 

dNjdp^ = " ^ (10) 

In pp collisions {nSi/ Sn)pp 1 due to the low string overlap probability. 

The STAR analysis of charged hadrons had presented the preliminary results for the 
percolation density parameter, ^ at RHIC for several collisions systems as a function of 
centrality[16]. Figure 2 shows ^ as function of the number of participant nucleons(iV)3art) 
in Au-|-Au collisions at ^/snn — 200 and 62.4 GeV. These experimental ^ values are used 
to get the bulk thermodynamic observables. 
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4. Determination of the Temperature 



The strong longitudinal chromo-electric fields produce Schwinger-Bialas [4,7,8] like ra- 
diation with a thermal spectrum, in analogy with the Hawking-Unruh radiation [17-22]. 
Both the Schwingcr-Bialas and Hawking-Unruh derivations lead to the same value of the 
maximum entropy temperature. Above the critical value of 5, the QGP in CSPM con- 
sists of massless constituents (gluons). The percolation density parameter ^ determines 
the cluster size distribution, the temperature T(^) and the transverse momentum in the 
collision [7]. The connection between ^ and the temperature T(^) involves the Schwinger 
mechanism (SM) for particle production. In CSPM the Schwinger distribution for mass- 
less particles is expressed in terms of 

dn/dp1 ~ e-'^P*/^' (11) 

with the average value of string tension {x^)- Gaussian fluctuations in the string tension 
around its mean value transforms SM into the thermal distribution [8] 

dn/dpi ~ e^-^'VW)^ (12) 
with {x^) = 7r(pj)i/-F(^). The temperature is given by 



rtt) = ^Ml (13) 

The string percolation density parameter ^ which characterizes the percolation clusters 
also determines the temperature of the system. In this way at ^c=l-2 the percolation 
phase transition at T(^c) models the thermal deconfincmcnt transition. In the determi- 
nation of temperature using Eq.(13) the value of F{^) is obtained using the experimental 
data [16]. We will adopt the point of view that the experimentally determined chemical 
freeze-out temperature is a good measure of the phase transition temperature, Tc [23]. 
The single string average transverse momentum {p1)\ is calculated at = 1-2 with the 
universal chemical freeze-out temperature of 167.7 ± 2.6 MeV [24]. This gives \J {pt)i 
= 207.2 ± 3.3 McV which is close to ~200 MeV used previously in the calculation of 
the percolation transition temperature [7]. Above =1-2 the size and density of the 
spanning cluster increases. We use the measured value of ^ = 2.88 to determine the 
temperature, before the expansion, Ti ~ 193.6±3.0 MeV of the quark gluon plasma in 
reasonable agreement with Tj = 221±19'**"* ± 19*^* from the enhanced direct photon 
experiment measured by the PHENIX Collaboration [2 5]. 

5. Bulk Thermodynamic Quantities 

Among the most important and fundamental problems in finite-temperature QCD are 
the calculation of the bulk properties of hot QCD matter and characterization of the 
nature of the QCD phase transition. The QGP according to CSPM is born in local 
thermal equilibrium because the temperature is determined at the string level. We use 
CSPM coupled to hydrodynamics to calculate energy density, pressure, entropy density 
and sound velocity. As mentioned earlier the strings interact strongly to form clusters 
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Fig. 3. The energy density 

and produce the pressure at the early stages of the colhsions, which is evident from the 
presence of elhptical flow in CSPM [14]. After the initial temperature T > Tc the CSPM 
perfect fluid may expand according to Bjorken boost invariant ID hydrodynamics [26] 



1 dT 

Tdr 

dT _ dT de 

dr de dr 



de 

d^ " 
. = (1 



-Ts/t 
de 



(14) 
(15) 
(16) 
(17) 
(18) 



where e is the energy density, s the entropy density, r the proper time, and Cg the sound 
velocity. Above the critical temperature only massless particles are present in CSPM. 
The initial energy density above T^ is given by [26] 



2 SjiT-nro 



(19) 



To evaluate ei we use the charged pion multiplicity dNc/dy at midrapidity and Sn values 
from STAR for 0-10% central Au-Au collisions with ^/JWn =200 GeV [27]. The factor 
3/2 in Eq.(19) accounts for the neutral pions. 

The dynamics of massless particle production has been studied in QE2 quantum elec- 
trodynamics. QE2 can be scaled from electrodynamics to quantum chromodynamics us- 
ing the ratio of the coupling constants [28]. The production time Tp^o for a boson (gluon) 
is [29] 



2A05h 

(mt) 



(20) 



In CSPM the total transverse energy is proportional to ^. From the measured value of ^ 
and e it is found that e is proportional to ^ for the range 1.2 < ^ < 2.88, = 0.788 ^ 
GeV/fm^ [16,27]. This relationship has been extrapolated to below ^ — 1.2 and above 
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Fig. 4. The speed of sound Fig. 5. The entropy density 

^ = 2.88 for the energy and entropy density calculations. Figure 3 shows e/T^ as obtained 
from CSPM along with the LQCD calculations [3] and the CSPM pressure Sp/T"^. 

For an ideal gas of massless constituents, the energy density and pressure are related 
by e = 3P. In LQCD the basic quantity is the interaction measure A = e-3p/T^, which 
is also known as the trace anomaly. In CSPM A reaches zero for T > 2Tc while in case 
of LQCD it is strongly interacting even for T > 4Tc [3] . 

The sound velocity requires the evaluation of s and dT/de, which can be expressed in 
terms of ^ and F{^). With g^/^ = F{^) one obtains; 



dT _ dT dq de_ 
de dq d^ de 

Then becomes: 



(21) 



= (-1/4)(1 + d) (^^^^ - 1 j (22) 

for ^ > ^c, an analytic function of £, for the equation of state of the QGP for T > Tc- 

Figure 4 shows the comparison of from CSPM and LQCD. The LQCD values were 
obtained using the EOS of 2+1 flavor QCD at finite temperature with physical strange 
quark mass and almost physical light quark masses [3]. At T/Tc—1 the CSPM and LQCD 
agree with the value of the physical hadron gas with resonance mass truncation M < 
2.5 GeV [30]. 

The entropy density s /T^ is obtained from energy density and speed of sound as shown 
in Fig. 5 along with the LQCD results. CSPM is in excellent agreement with the LQCD 
calculations in the phase transition region for T/Tc <1.5. 

6. Summary 

The percolation analysis of the color sources applied to STAR data at RHIC provides a 
compelling argument that the QGP is formed in central Au-Au collisions at y/s^lv = 200 
GeV. It also suggests that the QGP is produced in all soft high energy high multiplicity 
collisions when the string density exceeds the percolation transition. The results are also 
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in agreement with lattice QCD in the phase transition region, when the results are plotted 
with respect to T/Tf -SPA^ j^/j^lqcd^ j.-^^ ^^-^^^ =0.14 is in agreement with the 
physical resonance gas value at the critical temperature. Thus clustering and percolation 
can provide a conceptual basis for the QCD phase diagram which is more general than 
symmetry breaking [31]. 

7. Acknowledgement 

This research was supported by the Office of Nuclear Physics within the U.S. Depart- 
ment of Energy Office of Science under Grant No. DE-FG02-88ER40412. 



References 

[I] T. Cclik, F. Karsch, H. Satz, Phys. Lett. B97, 128 (1980). 
[2] M. B. Isichenko, Rev. Mod. Phys. 64, 961 (1992). 

[3] A. Bazavov et al., Phys. Rev. D80, 014504 (2009). 

[4] J. Schwinger, Phys. Rev. 82, 664 (1951). 

[5] M. A. Braun, C. Pajares, Eu. Phys. J. C16, 349 (2000). 

[6] M. A. Braun, F. del Moral, C. Pajares, Phys. Rev. C 65, 024907 (2002). 

[7] J. Dias dc Deus, C. Pajares, Phys. Lett. B642, 455 (2006). 

[8] A. Bialas, Phys. Lett. B466, 301 (1999). 

[9] H. Satz, Rep. Prog. Phys. 63, 1511 (2000). 

[10] J. Dias dc Deus, E. G. Ferreiro, C. Pajares, R. Ugoccioni, Phys. Lett. B581, 156 (2004). 

[II] J. Schaffner-Bielich, D. Kharzeev, L. McLerran, R. Venugopalan, NucL Phys. A705, 494 (2002). 
[12] J. Dias de Deus, E. G. Ferreiro, C. Pajares, R. Ugoccioni, Eur. Phys. J. C40, 229 (2005). 

[13] P. Brogueira, J. Dias de Deus, J. G. Milhano, NucL Phys. A. 832, 76 (2010). 
[14] L Bautista, L. Cunqueiro, J. Dias de Deus, C. Pajares, J. Phys. G37, 015103 (2010). 
[15] B. L Abelev et al., (STAR Collaboration), Phys. Rev. Lett. 103, 172301 (2009). 
[16] B. K. Srivastava, R. P. Scharenberg, T. Tarnowsky, (STAR Collaboration), Nukleonika 51, sl09 
(2006). 

[17] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975). 

[18] W. G. Unruh, Phys. Rev. D14, 870 (1976). 

[19] M. K. Parikh, F. Wilczek, Phys. Rev. Lett. 85, 5042 (2000). 

[20] D. Kharzeev, K. Tuchin, Nucl. Phys. A753, 316 (2005). 

[21] D. Kharzeev, E. Levin, K. Tuchin, Phys. Rev. C75, 044903 (2007). 

[22] P. Castorina, D. Kharzeev, H. Satz, Eur. Phys. J. C52, 187 (2007). 

[23] P. Braun-Munzinger, J. Stachel, Christof Wetterich, Phys. Lett. B596, 61 (2004). 

[24] F. Becattini, P. Castorina, A. Milov, H. Satz, Eur. Phys. J. C66, 377 (2010). 

[25] A. Adare et aL, (PHENIX Collaboration), Phys. Rev. Lett. 104, 132301 (2010). 

[26] J. D. Bjorken, Phys. Rev. D27, 140 (1983). 

[27] B. I. Abelev et al., (STAR Collaboration), Phys. Rev. C79, 34909 (2009). 

[28] C. Y. Wong, Introduction to high energy heavy ion collisions (World Scientific, 1994). 

[29] J. Schwinger, Phys. Rev. 128, 2425 (1962). 

[30] P. Castorina, J. Cleymans, D. E. Miller, H. Satz, arXiv:hep-ph/0906.2289vl. 
[31] H. Satz, Quantum field theory in extreme environments, Paris, April 2009. 



8 



